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Abstract: The well-known "Generalized Champagne Problem" on simultaneous stabilization of lin- 
ear systems is solved by using complex analysis [1, 11, 13, 18, 22] and Blondel's technique [5, 6, 8]. 
We give a complete answer to the open problem proposed by Patel et al. [20, 21], which automatically 
includes the solution to the original "Champagne Problem" [6, 8, 9, 11, 20, 21]. Based on the recent 
development in automated inequality-type theorem proving [31, 32, 33, 35, 36], a new stabilizing con- 
troller design method is established. Our numerical examples significantly improve the relevant results 
in the literature [11, 21 J. 
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1 Introduction 

Simultaneous stabilization of linear systems is a fundamental issue in system and control theory, and 
is of theoretical as well as practical significance [5, 15, 25] . The basic statement of the simultaneous 
stabilization problem of linear systems is as follows [5, 15, 23, 25]: 

Let pi,p2, ■ ■ ■ ,Pk be k scalar linear time-invariant systems. Under what condition does there 
exist a fixed controller c that is stabilizing for each pi(i — 1, ■ • • , k) ? 
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If k = 1, the above problem is reduced to the stabilization of a single system. There always exists 
a stabilizing controller for a single system provided no unstable pole-zero cancellations occur [5, 
15, 25]. Moreover, once a stabilizing controller of a single system is found, it is easy to parametrize 
the set of all stabilizing controllers of this system [5, 14, 15, 25]. This parametrization is known as 
Youla-Kucera parametrization discovered by Youla et al. [39, 40] and Kucera [16] respectively. 

When k = 2, by using the Youla-Kucera parametrization, it is possible to rephrase simultaneous 
stabilization of two systems into strong stabilization (stabilization with a stable controller) of a 
single system. This relationship was discovered for scaler systems by Saeks and Murray [23], 
and for multi- variable systems [5, 15, 25] by Vidyasagar and Viswanadham [26]. For the strong 
stabilization problem of a single system, Youla et al. [38] established an elegant criterion: a system 
is stabilizable by a stable controller if and only if it has an even number of real unstable poles 
between each pair of real unstable zeros! This remarkable and easily testable condition is called 
"parity interlacing property" [5, 15, 25, 38]. 

For k > 3, simultaneous stabilization problem is much more complicated than one expected 
[5, 6, 7, 8]. Vidyasagar and Viswanadham [26] stated that it is possible to transform simultane- 
ous stabilization of k systems to strong stabilization of corresponding k — 1 systems [5, 15, 25]. 
Blondcl et al. [8] proved that simultaneous stabilization of k systems is equivalent to bistable 
stabilization of associated k — 2 systems. Bistable stabilization means stabilization with a stable 
and inverse-stable controller. Such a controller is called bistable controller or unit controller [5]. 
Although many necessary or sufficient conditions for simultaneous stabilization of three or more 
systems were obtained in recent years (see [5] and references therein), easily testable necessary and 
sufficient conditions have not been found yet. Blondel and Gevers [6] showed that the simultane- 
ous stabilization of three systems is not rationally decidable [5, 6], i.e., it is not possible to find 
tractable necessary and sufficient conditions for simultaneous stabilization of three systems that 
involve only a combination of finite arithmetical operations (addition, substraction, multiplication 
and division), logical operations ('and' and 'or') and sign test operations (equal to, greater than, 
greater than or equal to, less than, less than or equal to) on the coefficients of the three systems! 

To illustrate the complexity of the simultaneous stabilization problem of three systems, Blondel 
et al. [6, 8, 9] proposed a specific simultaneous stabilization problem called "Champagne Problem" 
where the three systems are explicitly given. Patel [20] solved this problem by showing that there 
does not exist a stabilizing controller. Furthermore, a more general simultaneous stabilization 
problem, the "Generalized Champagne Problem", was proposed by Patel et al. [20, 21], which 
contains the "Champagne Problem" as a special case. 

The "Generalized Champagne Problem " is essentially concerned with determining the exact 
range of a specific parameter in simultaneous stabilization problem [20, 21]. Up to now, there are 
only some numerical estimates on this range in the literature [17, 21]. 

It should be pointed out that the "Champagne Problem" and the "Generalized Champagne 
Problem" are not of much engineering significance in themselves since they both deal with some 
specific systems. However, by studying these specific problems, we can develop new tools and 
new methods for general simultaneous stabilization problem. As a matter of fact, by studying 
some specific systems, Blondel et al. revealed the inherent relationship [5, 6, 8, 20] between the 
simultaneous stabilization problem and complex analysis theory [1, 11, 13, 18, 22], and derived the 
deep theoretic results of rational undccidability for three systems [5, 6]. 

In this paper, the well-known "Generalized Champagne Problem" on simultaneous stabilization 
of linear systems is solved by using complex analysis theory [1, 11, 13, 18, 22] and Blondel's 
technique [5, 6, 8]. We give a complete answer to the open problem proposed by Patel et al. 
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[20, 21], which automatically includes the solution to the original "Champagne Problem". Based 
on the recent advances in automated inequality-type theorem proving [31, 32, 33, 36, 35], a novel 
stabilizing controller design method is established. Our numerical examples significantly improve 
the relevant results in the literature [17, 21]. 

The paper is organized as follows. We state the generalized champagne problem and give a 
complete theoretical solution in Section 2. Section 3 deals with the controller design, and some 
numerical examples are provided. Finally, the conclusion is contained in Section 4. 

2 Solution to the Generalized Champagne Problem 

In this paper, all polynomials are of real coefficients. We denote by P the set of polynomials, 
P n the set of n-th order polynomials where n is an non-negative integer, H the set of Hurwitz 
stable polynomials (all roots lie within left half of the complex plane), H n the set of n-th order 
Hurwitz polynomials, S the set of Schur stable polynomials (all roots lie outside the unit circle) 1 
and S n the set of n-th order Schur polynomials. The variable of polynomials is s or z (we usually 
use s for Hurwitz polynomials and z for Schur polynomials). Let C denote the complex plane, 
D = {z G C : \z\ < 1} be the open unit disc and cl(D) = {z G C : \z\ < 1} be the closure of D, 
i.e., the closed unit disc. As usual, z stands for the complex conjugate of z G C, — C U {oo} 
the extended complex plane and C+oc = {s £ C : Rc(s) > 0} U {oo} the extended closed right half 
complex plane. 

Throughout this paper, plants and controllers are restricted to single-input single-output sys- 
tems that are described by linear, time-invariant, real rational transfer functions 2 . 

Given a system p, a controller c stabilizing p means that the four transfer functions pc(l + 
pc) _1 ,c(l + pc)^ 1 ,p(l + pc)~ x and (1 + pc)^ 1 are all stable, i.e., the denominator polynomials of 
the transfer functions are stable (which means Hurwitz stable in continuous-time case or Schur 
stable in discrete-time case). The degree of a controller is the maximum of the degrees of its 
denominator and numerator. 

viz) 

We denote by U (61(D)) = {c(z) = : x(z) G S,y(z) G S} the set of all bistable or unit 

controllers in discrete time case. Obviously, any element in U(cl(D)) and its inverse are analytic 
on cl(£>). 

The following problem is the well-known "Champagne Problem" [6, 8, 9, 17, 20, 21] of simul- 
taneous stabilization. It illustrates the difficulty in the simultaneous stabilization of three systems 
in general. It also shows that, simple problems do not always have simple answers. 

Champagne Problem [6, 8, 9, 17, 20, 21] Does there exist a controller that simultaneous 

1 Note that, some researchers use the concept of Schur polynomials for polynomials with all roots lying inside 
unit circle. Most papers on simultaneous stabilization adopt the definition used in this paper, which makes it easier 
to use the results in complex analysis. It is only a usage and doesn't affect the problem in essence. 

2 Here, the real rational transfer functions are not necessarily proper, that is to say, the degree of numerator may 
be greater than that of denominator. In the literature on simultaneous stabilization [5, 6, 8, 9], controllers do not 
have to be proper. It is mainly because that simultaneous stabilization problems can be transformed into the study 
on the properties of some analytic functions on the extend complex plane [1, 11, 13, 18, 22], In addition, Blondel 
[5] gave the conclusion that if k plants arc simultaneously stabilizable by a non-proper controller, then they are also 
simultaneously stabilizable by a proper controller. In fact, once a non-proper controller simultaneously stabilizes 
k systems, since the roots of a polynomial continuously depend on its coefficients [12, 37], a proper stabilizing 
controller can be obtained with a sufficiently small perturbation imposed on the denominator polynomial of the 
original controller. This technique will also be used in the following section on controller design. 
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stabilizes the following three plants: 

2 s - 1 , , (s - l) 2 , . „ n 

Pl(s) = T77TT' P2(s) = (9 S -8)(s + iy P3(s) = 0? 

Patcl [20] solved this problem in 1999 by showing that there does not exist a stabilizing con- 
troller. More generally, he considered the following problem: 

Generalized Champagne Problem [20, 21] What is the range of S for the existence of a 
controller that simultaneously stabilizes the following three plants: 

, \ 2£(s - 1) , . 2S(s - l) 2 , . 

where S is a real number. 

Obviously, "Generalized Champagne Problem" focuses on determining the range of 5 when there 
exists a simultaneously stabilizing controller while "Champagne Problem" asks whether S = — is 
in this range. 

Regarding "Generalized Champagne Problem" , as discussed in Patel [20] , there does not exist 

such a controller when 5 < — 3 . This conclusion certainly answers "Champagne Problem" and 

1.6 

Patel also stated that it remains an open problem [20] whether there exists such a controller when 
5 = — . Leizarowitz ct al. [17, 21] conjectured the minimum value of <5 for "Generalized Champagne 

Problem" is — = and they found a controller for this value of 5. But in 2002, Patel et 

2e 5.4366 

u { s ) 

al. [21] showed that there exists a stabilizing controller c = — — where x(s) £ H 9 ,y(s) G P 9 

x(s) 

when S = - — — — — Note that the degree of the controller therein equals 9. Hence, the 

6.719367588932806 6 H 

conjecture mentioned in [17] is invalid. Thus for "Generalized Champagne Problem", it is still an 

1 1 

open problem whether there exists a stabilizing controller when 5 £ \ — , ). 

6 L 16 6. 719367588932806 ; 

For "Generalized Champagne Problem" with S — — , since the roots of a polynomial continu- 

16 

ously depend on its coefficients [12, 37], it is easy to know that there does not exist a controller that 
simultaneously stabilizes the three plants. Thus we have given an answer to the problem proposed 
in Fatal [20]. 

It should be pointed out that the above results presented in Patel et al. [21] and Leizarowitz 
et al. [17] are both based on numerical analysis. The degrees of the controllers provided are high 
and the conclusions are not theoretically complete. 

The following theorem is the main result of this paper which gives a completely theoretical 
solution to "Generalized Champagne Problem" . 

Theorem 1 The necessary and sufficient condition for the existence of a controller that 
simultaneously stabilizes the following three plants: 

2S(s-l) 25(s - l) 2 

Pl(s) = ^TT-' P2{s)= ((i + 5) s -(i-6))( s + iy P3(s) = (1) 

where 5 is a real number, is S = or \S\ > — . 
16 

3 In fact, Patcl [20] admitted <5 > 0. Otherwise this condition should be < |<5| < — . Indeed, from the result 

1 ^ 

in the sequel, when <5 = Oor<5<— — , there does exist a simultaneously stabilizing controller for "Generalized 

Champagne Problem". 
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To prove the above theorem, the following two well-known results in complex analysis will be 
used. ^ 

Lemma 1 (jq — Theorem) [5, 13, 18, 20] 4 Any analytic function f on D that satisfies: 

1) f(z) =0,zeD<»z = 0, 

2) /'(0) = 1, 

contains an open ball centered at z = with radius — in its range on D but not always a larger 
ball. 

Lemma 2 (Runge's Theorem) [11, 22] If £1 is an open subset of C such that \ £1 is 
connected then for each analytic function f onVL there is a sequence of polynomials {q n } such that 
q n uniformly converges to f on any compact subset of Cl. 

We are now in the position to give a proof of Theorem 1 . 

1 + z 

Proof of Theorem 1 The well known linear fractional transformation s = is a one-to- 

1-2 

S — 1 

one mapping of cl(D) onto C+oo and its inverse is z — . By this transformation, "Generalized 

s + 1 . . 

Champagne Problem" is equivalent to the simultaneous stabilization problem of the following three 
plants in discrete time case: 

25z 2 

Pid(z) = 2Sz, p 2d {z) = — —z, P3d(z) = 0, (2) 

2 + 

where 6 is a real number. 

The proof contains the following three steps. 

Step 1 Firstly, we prove that the three plants defined by Eq. (2) are not simultaneously 

stabilizable in discrete time case when < \S\ < — . 

16 ^ 

Assumed by contradiction that, for < \6\ < — , there is a stabilizing controller c(z) for 

the three plants pid(z) , Pid(z) , pzd(z) . Then, since c(z) stabilizes pzd(z) = 0, c(z) must be stable. 

25z 2 

Moreover, c(z) also stabilizes pid(z) — 2Sz and pid{z) = — - and we have the following equations: 

2Szc(z) + l = Ul (z) e U(cl(D)) (3) 
2Sz 2 c(z) + z + S = u 2 (z) E U (61(D)) (4) 

Set 

f(z) = 25z 2 c(z) + z = z(26zc(z) + 1). 
Obviously, f(z) is analytic on D and f(0) = as well as /'(0) = 1. Also, we have 

f(z) = z Ul (z) (5) 
f(z) + 5 = u 2 (z) (6) 

The Eq. (5) means that in D, f(z) is equal to if and only if 2 is while the Eq. (6) implies 
that —5 is not in the range of f(z) in D. However, f(z) satisfies the hypothesis of — — Theorem, 

but < \5\ < — . This is impossible, a contradiction is obtained and Step 1 is proved. 

Step 2 When (5 = 0, the three plants defined by Eq.(2) degenerate to a plant p\d(z) = 
Pid(z) = Pzd(z) = 0. Apparently, c(z) meets the requirement if c(z) is stable. 

When \6\ = — , since the roots of a polynomial continuously depend on its coefficients [12, 37], 
16 

it is easy to know that there does not exist the desired controller, namely, the three plants defined 
by Eq. (2) are not simultaneously stabilzable in this case. 

4 There is a misprint in [20] for the reference of this lemma where Condition 2 is misprinted as /(0) = 1. 
5 This lemma is a special case of Runge's Theorem, i.e., a corollary of the original Runge's Theorem [11, 22]. 
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Step 3 We prove that the three plants defined by Eq.(2) are simultaneously stabilizable in 

discrete time case when \S\ > —. 
1 

If | S | > — j from Lemma 1, there exists an analytic function f(z) on D such that f(z) = 
4i> z = and /'(0) = 1 but — ^ f(D). Moreover, f(z) can be a real analytic function, i.e., 
f(z)=W)- 6 



Set 



9{z) = - (7) 
According to the properties of f(z), it is obvious that g(z) satisfies the following properties: 
!) 9(z) = 9{z), 

2) g(z) ia analytic on \z\ < 16\S\, 

3) g(z) = z = 0, and g'(0) = 1, 

4) 5(2) leaves out —S on z < 16\S\. 
Let 

M = mf |<?(z) + <5|. (8) 

z£Cl{D) 

and 

From the above properties of g(z), we have \i > 0, is analytic on \z\ < 16|£|, and 

1 + 26z/i(z) ^ 0, zecl(D). (10) 

By Runge's Theorem, there exists a real polynomial q(z) such that 

\h(z)-q(z)\<-±(±) 2 , zec\(D). (11) 

Since h(z) satisfies (10), q(z) also satisfies 

1 + 2Szq(z) ^ 0, zecl(D), (12) 

i.e., 

l + 26zq(z) = u 1 (z) e U(cW). (13) 
By the definition of h(z), we have 

\g(z)- z-2Sz 2 q(z)\ < fi, z e cl(D). (14) 

Set 

p(z) =z + 2£z 2 <?(z). (15) 
Obviously, p(0) = 0,p'(0) = 1. By the definition of /z, 

p(z) 96 —5, z G cl(D). (16) 

Hence, we have 

5 + z + 2Sz 2 q(z) =u 2 (z) e U(cW). (17) 



6 



Special constructions of such functions are presented in [13, 18], for instance, 

+ z 2n 

T 

n=l 



«*>=*II<it5^> 8 
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Let c(z) = q(z), by (13) and (17), c(z) is the desired controller. That completes the proof of 
Step 3. 

Theorem 1 is now proved by combining Steps 1-3. 

Remark 1 Theorem 1 contains a complete theoretical answer to the 'Generalized Champagne 
Problem " [20, 21], automatically including the solution to the "Champagne Problem" [6, 8, 9, 20]. 
The main result in Patel [20] is also contained in Step 1 of the proof of Theorem 1. Apparently, 
the proof in this paper is much more concise and straightforward than that in Patel [20]. 

Remark 2 As a byproduct, Step 2 of the proof of Theorem 1 contains an answer to an open 
problem proposed in Patel [20]. 

By a similar argument as in the proof of Theorem 1, we have 

Theorem 2 The necessary and sufficient condition for the existence of a controller that 
simultaneously stabilizes the following three plants in continuous time case: 

^ = 7TI< »M= «i-s)*-v!-{i + s) > P3(s) = °< (18) 

where S is a real number, or, equivalently, simultaneously stabilizes the following three plants in 
discrete time case: 

z 2 

pid(z) = z, p 2 d{z) = 7, p 3 d(z) = 0, (19) 

z — o 

where S is a real number, is S — or \S\ > — . 

1 1 16 

Remark 3 Some partial results of Theorem 2 are available in [5, 7, 8]. In [5[, a sufficient 
condition was given on the nonexistence of a simultaneously stabilizing controller for the three 
plants (in continuous time case) in Theorem 2. 

Remark 4 It should be pointed out that Theorem 2 is equivalent to Theorem 1 in some sense. 
In fact, a stabilizing controller multiplied by a non-zero constant is also a stabilizing controller. 
When 5 7^ 0, by multiplying the three plants in (19) by —25, and letting 5\ = —8, it is obvious that 
the three plants obtained are the same as the three plants in Theorem 1. 

Remark 5 Although a complete theoretical solution to the "Generalized Champagne Problem" 
is provided in Theorem 1, it can be seen from the proof of Theorem 1 that it is still difficult in 
practice to construct the simultaneously stabilizing controller. This difficulty can also be seen by the 
controller design examples in [17, 21]. A novel stabilizing controller design method is established 
in the next section, and the numerical examples obtained significantly improve the relevant results 
in [17, 21]. 



3 Controller Design 

Although the "Generalized Champagne Problem" has been theoretically resolved, it is still diffi- 
cult in practice to construct the simultaneously stabilizing controller. This fact can also be seen 
in the specific controller design examples in [17, 21]. In particular, how can one construct the 

simultaneously stabilizing controller with minimal degree for d > — ? And when the degree of 

16 

the controller is fixed, how to determine the range of 67 These issues are of both theoretical and 
practical significance. 

When the degree of the controller is fixed, the controller design problem of simultaneous stabi- 
lization can be transformed in essence to the problem of how to solve a set of algebraic inequalities. 
Early in 1950s Tarski published the well-known work [24] on the decidability [2, 24, 27, 28, 30, 37] 
of this kind of problems. Tarski's decision algorithm is of theoretical significance only, since it can 
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not be used to verify any non-trivial algebraic or geometric propositions in practice [27, 28, 30, 37] 
due to its very high computational complexity. The Cylindrical Algebraic Decomposition (CAD) 
[3, 4, 10] algorithm proposed by Collins and subsequently improved by him and his collaborators 
is the first practical decision algorithm and can be used to verify non-trivial algebraic or geometric 
propositions on computer. Although, as a generic program for automated theorem proving, its 
computational complexity was still very high, the CAD and its improved variations have become 
one of the major tools for solving this kind of problems. 

Wu [27] proposed in 1977 a new decision procedure for proving geometry theorems. As an 
important progress in automated theorem proving [28, 30, 37], Wu's method is very efficient for 
mechanically proving elementary geometry theorems of equality type. The success of Wu's method 
has inspired in the world the research of algebraic approach to automated theorem proving [28, 
30, 37]. However, automated inequality proving has been a difficult topic in the area of automated 
reasoning for many years since the relevant algorithms depend on real algebra and real geometry. 
In 1996, Yang and his colleagues [34, 37] introduced a powerful tool, the Complete Discrimination 
Systems (CDS) of Polynomials, for automated reasoning in real algebra. By means of CDS, many 
inequality-type theorems from various applications have been proved or discovered. 

In recent years, by combining discriminant sequences of polynomials [34, 37] with Wu's method 
[27, 28, 30] as well as Partial Cylindrical Algebraic Decomposition [3, 4, 10], Yang et al. proposed 
some algorithms which are able to discover new inequalities. These algorithms arc complete for 
an extensive class of problems involving inequalities and are applicable to the controller design 
in simultaneous stabilization. Based on these algorithms, a generic program called "Discoverer" 
[34, 35] and a generic program called "Bottema" [31, 32, 36] were implemented as Maple packages. 
The following controller examples are all obtained automatically by calling "Discoverer" [34, 35] 
or "Bottema" [31, 32, 36]. 

The following sufficient condition for checking the Hurwitz stability of a polynomial may be 
used for improving computational efficiency when calling the above programs. 

Lemma 3 [19] Suppose f(s) = a Q + a x s + • • • + a n s n (a^ > 0, i = 0, 1, 2, ■ • • , n, n > 3), f(s) 
is Hurwitz stable if 

ai-\o,i+2 < 0.4655aia i+ i, (i = 1, 2, • • • , n - 2). (20) 

In addition, to get a proper controller, we introduce the following theorem on the disturbance 
to the Hurwitz polynomials. 

Theorem 3 Given a polynomial f(s) = a + a\s + • • • + a m s m (<Zj > 0, i = 0, 1, 2, • ■ ■ , to) G 
H m , then for any integer n > m and sufficiently small e > 0, there exist e i; < Si < e (i = 
to+1,to+2, • • • , n), such that g(s) = a^+ais^ ^a m s m +e m+1 s m+1 +e m+2 s m+2 -\ \-e n s n E H n . 

Proof Suppose f(s) = a + ais H — • + a m s m (ai > 0, i = 0, l, 2, ■ • • , m) € H m , we only need 
to prove that g(s) = a + a^s + • • • + a m s m + es m+1 6 H m+1 for sufficiently small e > because 
we can treat g(s) with the same manner iteratively. 

Since f(s) = a + ais + • • • + a m s m G H m , it is well known that Hurwitz principal minors 
A{ >0(i = l,2,.--,m) [12]. 

By computing Hurwitz principal minors of g(s) = oq + a\s + ■ ■ ■ + a m s m + es m+1 , it is 
easy to verify that — a m , A? = a m A{_ 1 + e§i(e, ao, a±, ■ ■ ■ , a m ) (i — 2, ■ • • , m, m + 1), where 
5i(s, ao, oi, • • • , a m ) (i — 2, • • ■ , to, m + 1) are polynomials in (e, a , a\, ■ ■ ■ , a m ). Obviously, when 
e > is sufficiently small, we have A^ > (i = 1, 2, • • • , to, to + 1). This completes the proof. 

In what follows, we will provide some controller examples for "Generalized Champagne Prob- 
lem" of simultaneous stabilization by applying "Discoverer" [34, 35] or "Bottema" [31, 32, 36]. 
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Without loss of generality, we suppose the denominator polynomials of the controllers are monic. 
In addition, for convenience, we only consider the case 5 > 0. 

Example 1 For the three plants defined in "Generalized Champagne Problem", when the 
degree of controllers is restricted to 0, i.e., the controllers are constants, invoking "Discoverer" or 
"Bottema", we obtain the following results. There exist desired controllers for S > —, whereas there 

113 

do not exist requested controllers for 8 < —. For a given 8 > —, for example, 8 = —, c(s) = yo is 

a desired controller if and only if — < y < —. 

Remark 6 Example 1 looks simple, but it shows that: 1 ) When the degree of the controller 
is restricted, the explicit bound of 8 can be obtained by using the packages developed by Yang et al; 
2) Moreover, when 8 is fixed, the ranges of the parameters of controllers can be obtained. Due to 
the completeness of these algorithms, the conditions obtained are both necessary and sufficient. 

Example 2 For the three plants defined in "Generalized Champagne Problem", when the 
numerator polynomials of controllers are restricted toO — th order polynomials and the denominator 
polynomials of controllers to 1st order polynomials, i.e., the controllers have the form c(s) = 
— — — , invoking "Discoverer" or "Bottema", we have the following results as in the above example. 

S + Xo 

There exist desired controllers for 8 > —, whereas there do not exist requested controllers for 8 < —. 

13 1 

For a qiven 5 > —, say, 8 = —, c(s) = is a desired controller. 

y 2' y ' 4' w s + 3 

Example 3 For the three plants defined in "Generalized Champagne Problem", when the 

numerator polynomials of controllers are restricted to 1st order polynomials and the denominator 

polynomials of controllers to — th order polynomials, i.e., the controllers have the form c(s) = 

VI s + yo, invoking "Discoverer" or "Bottema", we have the following results. There exist desired 

1 1 
controllers for 8 > -, whereas there do not exist requested controllers for 8 < -. For a given 

S> \> sa V S = \> ifiyi'Vo) e {(-^, (jj, jj), (^, jj)}, c(s) = yis + yo is a requested controller. 

Moreover, to get a proper controller, by the continuous dependance for roots of polynomials on their 

151 3 3 4 3 

coefficients [12, 37] and Theorem 3, if {x 1 ,y 1 ,y ) E {(e, -77,77^7:), (e, ■=, t), (e, 7, 7)} and e > is 

lu 1UU 5 5 5 5 

sufficiently small, e.g., e — — , c(s) = — — is a desired proper controller. 

10 es + 1 

Example 4 For the three plants defined in "Generalized Champagne Problem", when the 
numerator polynomials of controllers are restricted to 1st order polynomials and the denominator 
polynomials of controllers to 1st order polynomials, i.e., the controllers have the form c(s) = 

— — ^° , invoking "Discoverer" or "Bottema", we have that there exist desired controllers for 
s + x 

1 111 

8 > — whereas there do not exist requested controllers for 8 < —. For a given 8 > -, say, 8 = -, 

if(xo,yi,Vo) G {(2, ^, ^1), (3, ^, ^J), (6, 8, ^)}, c(s) = is a desired controller. 

Remark 7 The controllers obtained in the above examples are the sample points picked out 
from the cells of some appropriate decomposition of the parametric space which satisfy the require- 
ments of simultaneous stabilization. As mentioned in Example 1, when the degrees of controllers are 
restricted, the explicit bound of 8 and the ranges of the parameters of the controllers can be obtained 
by further analysis on the distribution of the roots of the polynomials output by "Discoverer" . 

Example 5 For the three plants defined in "Generalized Champagne Problem", when the 
numerator polynomials of controllers are restricted to 2nd order polynomials and the denominator 
polynomials of controllers to — th order polynomials, i.e., the controllers have the form c(s) = 
J/2S 2 + Vis + yo, invoking "Discoverer", we have that there do not exist desired controllers for 
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S = — . When 8 = — , using the sufficient condition proposed in Lemma 3 for Hurwitz stability, 
6 59 

191 39001 

we have that there exist such kind of controllers. For example, if (y2,Vi ,Vn) £ I ( , , 

j w*i>iw i v 100 , 10000 ' 

2450001 191 390019 2450003 97 39001 245001 97 19501 245001 19501 39003 

1000000^' ^100' 100000' 1000000 ^50' 10000' 100000 ^50' 5000 ' 100000^' U0000' 10000' 
4900299,, , ,.,„,, „ , 

2000000 C = ^ + 2/i- s + Vo ls a desired controller. Moreover, to get a proper controller, by 

the continuous dependance for roots of polynomials on their coefficients [12, 37] and Theorem 3, if 

191 39001 2450001 191 390019 2450003 97 39001' 

{x 2 ,x u y2,yi,yo) G {(£, e, —, 1000000 )> ( e , e > ^ob' TooOOO' 1000000 (£ ' £ ' 50' 10000' 

245001, , 97 19501 245001, , 19501 39003 4900299, , 

), (£, £, — , , ), (£, £, , , ) \ and e > is sufficiently 

100000 ; ' v ' '50' 5000' 100000 ; ' v ' ' 10000' 10000' 2000000" JJ y 

small, e.g., e — then c(s) = ^ 2 „ — — is a desired proper controller. 

y 10000000 w es 2 + £.s + l 

Remark 8 The value of S obtained in Example 5 is an improvement over the bound conjec- 
tured in [17]. Although an improvement over this bound was also made by Patel et al. [21], the 
degree of the controller provided by them equals 9 whereas the controller presented in Example 5 
is of degree 2 which is much more lower. In addition, we can see from the above examples that 
the improvement over the bound of S conjectured in [17] could not be achieved by controllers with 
degrees less than 2. 

Example 6 For the three plants defined in "Generalized Champagne Problem", suppose 
the numerator polynomials of controllers are restricted to 3rd order polynomials and the de- 
nominator polynomials of controllers to — th order polynomials, i.e., the controllers have the 
form c(s) = y 3 s 3 + y 2 s 2 + y\S + yo- When 6 = —, using the sufficient condition proposed in 

Lemma 3 for Hurwitz stability, we have that there exist such kind of controllers. For example, if 

1037 30077 50001 300001 26 376 50001 300001 113 378 50001 
(2/3,2/2,2/i,2/o) G { ( 1000 > Y0000' 10000' 100000'*' ^25' 125' 10000' 100000^' ^100' 125' 10000' 
300001 57029 121 50001 300001 11407 121 50001 300001 114113 1513 50001 

100000^' *■ 50000' lb~' 10000' 100000'*' 10000' lb"' 10000' 100000 "' *■ 100000' 500 ' 10000' 
300001,, , , 

100000 = V 3 V 2 y iS ^y° ls a desired controller. Moreover, to get a proper controller, 

by the continuous dependance for roots of polynomials on their coefficients [12, 37] and Theorem 3, 
... . ,. 1037 30077 50001 300001, , 26 376 50001 

if (^3,^2,^,2/3,2/2,2/1,2/0) G {(£1, £, £, — , — , — , ^- Q ), £, £, -, — , — , 

300001 113 378 50001 300001 57029 121 50001 300001 

looooo'*' ^ £l ' £ ' £ ' Too' 125' ioooo' iooooo^' ^ £l ' £ ' £ ' 50000' lb"' ioooo' looooo'*' ^ £l ' £ ' £ ' 

11407 121 50001 300001, , 114113 1513 50001 300001,, , , „, . , 

, , , ), (£i , £, £, , , , ) } and £ > £i for sufficiently 

10000' 40 ' 10000' 100000^' v ' ' ' 100000' 500 ' 10000' 100000 ;/ 2 

small £i > and e > 0, e.g., e — 10~ 7 and e\ = 10~ 15 , c(s) = ^ 3S „ 2/2 — _i yo - g a ^ es j re ^ 

eis-s + es z + es + 1 

proper controller. 

Remark 9 The value of 6 appeared in Example 6 is an improvement over the minimal bound 
proposed in [21]. The degree of the controller presented in [21] equals 9 whereas the controller given 
in Example 6 is of degree 3 which is much more lower. 

Remark 10 The above algorithms are also complete for high-order controller design. From 
our computational experiments, other phenomena can be observed. For instance, the improvement 
on the bound of 5 in "Generalized Champagne Problem" mainly depends on the increase on the 
order of numerator polynomial of the stabilizing controller. These problems deserve further research 
and are omitted here for succinctness. 

Remark 11 The generic programs "Discoverer" [33, 35] and "Bottema" [31, 32, 36] are 
powerful tools in practice for automated inequalities proving and can be applied in various fields. 
In this paper, only some basic functions of those packages are employed to determine the ranges of 
parameters and to design simultaneously stabilizing controllers. In previous examples, it is demon- 
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strated how powerful these packages are. Indeed, according to the Finitcncss Theorem proposed by 
W.T. Wu [29] on global-optimization, global- optimization problems can be theoretically solved by 
'Discoverer" [33, 35] and "Bottema" [31, 32, 36] where the objectives are polynomials and con- 
straints are also polynomial equalities or inequalities. The potential applications of the packages 
are considerable. 

Remark 12 It should be pointed out that, although "Discoverer" [33, 35] and Bottema [31, 
32, 36] are powerful, the computational complexity increases very quickly with the dimension, i.e., 
the number of parameters. It is a problem deserving further study how to promote the efficiencies 
of these algorithms for dealing with the large-scale engineering optimization problems. Combining 
symbolic computation with numerical calculation as well as large-scale parallel numerical processing 
may be an effective way [30]. 

4 Conclusion 

In this paper, the well-known "Generalized Champagne Problem" for simultaneous stabilization 
of linear systems has been resolved by using complex analysis [1, 11, 13, 18, 22] and Blondel's 
technique [5, 6, 8]. We gave a complete answer to the open problem proposed in Patel et al. 
[20, 21], which automatically includes the solution to the original "Champagne Problem" [6, 8, 
9, 17, 20, 21]. Based on the recent developments in automated inequality-type theorem proving 
[31, 32, 33, 36, 35], a novel stabilizing controller design method has been established. Our numerical 
examples significantly improved the relevant results in the literature [17, 21]. 
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